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ABSTRACT 

A detailed analysis of the data from a high sampling rate, multi-month reverberation mapping 
campaign, undertaken primarily at MDM Observatory with supporting observations from telescopes 
around the world, reveals that the II/3 emission region within the broad line regions (BLRs) of several 
nearby AGNs exhibit a variety of kinematic behaviors. While the primary goal of this campaign was 
to obtain either new or improved II/3 reverberation lag measurements for several relatively low lumi- 
nosity AGNs, we were also able to unambiguously reconstruct velocity-resolved reverberation signals 
from a subset of our targets. Through high cadence spectroscopic monitoring of the optical contin- 
uum and broad H/3 emission line variations observed in the nuclear regions of NGC3227, NGC3516, 
and NGC 5548, we clearly see evidence for outflowing, infalling, and virialized BLR gas motions, 
respectively. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert 
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1. INTRODUCTION 

Reverberatio n mapping (IBlandford fc McKed [l98l 
iPetersonI [l993l ) is a technique applied to spectroscopic 
observations of type 1 AGNs to infer properties of the 
broad line-emitting region (BLR) through characteri- 
zations of time delays between continuum and broad 
emission-line flux variations. This method has become 
extremely useful and quite successful in its current ap- 
plication of directly measuring BLR radii and black hole 
masses (Mbh)- However, the fundamental objective of 
reverberation mapping, as its name implies, is to re- 
construct or map the emissivity and velocity distribu- 
tion of the BLR line-emitting gas as a function of po- 
sition as it 'reverberates' in response to the flux varia- 
tions of the ionizing continuum. Because the position of 
the gas can be inferred through the emission-line time 
delay, r (i?BLR = cr), the res ulting reconstructi on is 
called a velocity-delay map (see lHorne et al.ll2004f ) and 
is the best means, with current technology, to obtain di- 
rect knowledge about the geometry and kinematics of the 
BLR. 

Past attempts at producing velocity-delay maps (e.g.. 
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Done fc Krolig[l99l: lUlrich fc Honii[l99l : iKollatschnvl 
2003f) have not yielded completely satisfactory results, 



primarily on account of limitations in temporal sam- 
pling, with inadequate time resolution or campaign du- 
ration or both. Even given these limitations, previ- 
ous studies have successfully measured time delays and 
black hole masses in over 40 typ e 1 AGNs (see e.g., 
iPeterson et"all 12004 [Bentz et al.l I2OO90 . In addition, 
basic investigations into time delay differences between 
multiple emission lines have shown that the BLR is virial- 
izcd ac ross broad line- emit ting regions of different species 
(e.g. JPeterson et ai]f2004 and references therein). Other 
studies of the velocity dependence of the lag across a sin- 
gle emission line region have shown suggestive evidence 
that the BLR commonly contains a radial inflow compo- 
nent in addition to circular motions (e.g., GaskcU 1988; 
Crensh aw fc Blackwelll [iQQOt iKoratkar fc Gaskelll 1199X1 : 
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The experiences from earlier reverberation programs 
have led to more recent campaigns that address the main 
observational obstacles encountered in the past — rela- 
tive sampling rate, campaign duration, and data qual- 
ity. Consequently, there has been consistent success in 
measuring BLR radii in new targets and remeasuring 
radii that now supersede previous, lower-precision or am- 
biguo iis measurements due to inadequate tirne-sampling 



„ 

(e.g. iDennev et al.l [20061: iBentz et al.l 120061: iGrier et all 



120081 : Dcnncy et all l200a hereafter D09a). Further- 



more, these campaigns enable statistically significant de- 
tections of reverberation signals at higher velocity res- 
olutions than have previously been found due mainly 
to the data restri ctions discussed above. For example, 
IBentz et al.l (|2008f ) revealed a clear signal of radial inflow 
within the H/3 em ission region of the BLR in Arp 151, 
and Bcntz 'et al.l (j20d9) show further evidence for dis- 
tinct kinematic behavior in at least one other AGN from 
the same program (the Lick AGN Monitoring Project). 
Here we report on results of our own recent high time- 
resolution monitoring program undertaken at five ob- 
servatories. The first results of this program have al- 
ready been published (D09a), and additional results are 
in preparation (K. D. Denney et al. 2009, in prepara- 
tion, hereafter D09b). In this Letter, we present the 
most clear velocity-resolved reverberation signals recov- 
ered from this campaign. Our results highlight three 
distinctly different BLR kinematic signatures — inflow, 
outflow, and virialized motions — for three separate tar- 
gets — NGC3227, NGC3516, and NGC5548, respec- 
tively. These results reveal kinematic diversity in rever- 
beration signals and underscore the importance of higher 
time resolution spectral monitoring. This represents an 
important step between measuring average BLR response 
times and black hole masses to realizing the full potential 
of this technique through the recovery of a velocity-delay 
map. 

2. OBSERVATIONS AND DATA ANALYSIS 

Spectra of the nuclear regions of NGC 3227, NGC 3516, 
and NGC 5548 were obtained from a combination of 
the 1.3 m telescope at MDM Observatory, the 2.6 m 
Shajn telescope of the Crimean Astrophysical Observa- 
tory (CrAO), and the Plaskett 1.8 m telescope at Domin- 
ion Astrophysical Observatory (DAO). Spectroscopic ob- 
servations targeted the H/3 A4861 and [O ill] AA4959, 5007 
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Fig. 1.— Mean and rms spectra of NGC 3227 (left), NGC 3516 
(middle), and NGC 5548 (right) from MDM observations. The 
solid line shows the mean and rms spectrum formed after removal 
of the [O in] AA4959, 5007 narrow emission lines, and the dotted 
lines represent the spectra prior to this subtraction (note the small 
residuals in the rms spectra). 

emission line region of the optical spectrum. The top 
panels of Figure [T] show the mean spectrum of each of 
the three targets based on the MDM observations, and 
the bottom panels show only the variable emission in the 
form of an rms spectrum for each object, respectively. 
Emission line light curves were made from the integrated 
H/3 flux measured above a linearly interpolated contin- 
uum fit to the MDM, CrAO, and DAO spectra. 

In addition to spectral observations, we obtained sup- 
plemental y-band photometry from the 2.0-m Multicolor 
Active Galactic NUclei Monitoring (MAGNUM ) tele- 
scope at the Haleakala O bservatories in Hawaii (|Yoshiil 
120021: lYoshii et al] 120031) . the 70-cm telescope of the 
CrAO, and the 0.4-m telescope of the University of Ne- 
braska (UNebr.). Continuum light curves were created 
with observations from each U-band photometric data 
set and the average continuum flux density near rest 
frame ~5100 A in each spectrum of the spectroscopic 
data sets. The reader is referred to D09a and D09b for 
details describing campaign observing setups and data 
reduction, flux calibration of the spectra, and intercali- 
bration of the data sets to form the single set of optical 
continuum and H/3 emission line light curves shown for 
each object in Figure [2l 

Table [T] displays basic statistical parameters describing 
the final light curves shown in Figure O Column 1 gives 
the object, and Column 2 lists the spectral feature rep- 
resented by each light curve. The number of data points 
in each light curve is shown in Column 3, with the me- 
dian sampling interval between these data points given 
in Column 4. Column 5 shows the mean fractional er- 
ror in the fluxes of each time series. Column 6 gives the 
excess variance, calculated as 
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where is the variance of the observed fluxes, 5^ is their 
mean square uncertainty, and (/) is the mea n of the ob- 
served fluxes (jRodriguez-Pascual et al.l[l997l ) . Column 7 
is the ratio of the maximum to minimum flux in the light 
curves, and Column 8 gives the adopted mean H/3 time 
lag and uncertainties determined through the primary 
time series analysis which utilized the full H/3 line profile 
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TABLE 1 

Light Curve Statistics and Mean H/3 Lags 



Object 
(1) 


Time 
Series 
(2) 


N 
(3) 


median 

(days) 
(4) 


Frac. Err 

(5) 


Mean 
(6) 


_Rmax 
(7) 


Tccnt (days) 
(8) 


NGC 3227 
NGC 3516 
NGC 5548 


5100 A 

H/3 
5100 A 

H/3 
5100 A 


171 
75 
198 
93 
182 
108 


0.45 
1.00 
0.54 
1.00 
0.56 
1.00 


0.03 
0.03 
0.06 
0.04 
0.03 
0.09 


0.10 
0.08 
0.28 
0.15 
0.11 
0.26 


1.9±0.1 
1.5±0.1 
5.9 ±1.5 
1.9 ±0.2 
1.7±0.1 
3.7±0.5 


3.8 ±0.8 
ll.7ll;» 

12.41^3:^ 
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Fig. 2. — Final light curves, after merging all data sets, of the 
5100A continuum (top panels) and broad H/3 emission line flux 
(bottom panels) in units of 10~^^ erg s~^ cm~^ A~^ and 10~^^ erg 
s-l cm-2, respectively, for NGC 3227 (top), NGC 3516(middle), 
and NGC 5548 (bottom). Except for the continuum flux scale of 
NGC 3227, which is arbitrary because a linear fit to long-term sec- 
ular trends has been divided out, the flux scales reflect all relative 
and absolute flux calibrations described by D09b. 

and is described in detail for each object by D09b. 

3. VELOCITY-RESOLVED TIME SERIES ANALYSIS 

The lag measurements between the continuum and H/3 
emission listed in Table [1] represent the average time de- 
lay across the BLR, measured from the centroid o f the 
cross correlation function fsee iPeterson et aLlll998l and 
references therein) produced during the time series anal- 
ysis of the continuum and full H/3 line profile light curves 
shown in Figure (see D09a for details). Here, we focus 
on the velocity-resolved time series analysis we performed 
on each target to investigate the potential for recovering 
velocity-dependent time delays across the H/3 emission 
line in order to infer the kinematic structure of the line 
emitting gas. 

We divided the H/? emission line into eight velocity- 



space bins, whose boundaries were determined by the 
division of the rms spectrum of each object into eight bins 
of equal flux, as depicted in the top panels of Figure [3l 
Compared to the line boundaries used for the full profile 
analysis leading to the light curves shown in Figure [21 
those used for this analysis were slightly narrowed in the 
cases of NGC 3516 and NGC 3227 in order to include 
only the most variable portions of the line profile, and 
boundaries were broadened for NGC 5548 because the 
rms spectrum shows variability in the red H/3 wing that 
extends beneath the [O ill] A4959 narrow emission line. 

Light curves were created from measurements of the 
integrated H/3 flux in each bin and then cross corre- 
lated with the continuum light curves shown in Figure 
[21 The bottom panels of Figure [H show the lag measure- 
ments for each of these bins. Error bars in the veloc- 
ity direction represent the bin width. The evidence for 
a velocity-stratified BLR response to continuum varia- 
tions is clear in all three cases. Interestingly, each case 
demonstrates a different kinematic signature: (1) out- 
flow is indicated in NGC 3227, given the generally longer 
lags from red-shifted BLR gas compared to the gas on 
the blue-shifted side of the line, (2) NGC 3516 shows 
the opposite signature, with the blue side of the line 
lagging the red side — an indication that there is an 
infall component to the gas i n this region, simi lar to 
that observed in Arp 151 bv iBentz et al.l ()2008f ). and 
(3) NGC 5548 shows no radial gas motions, with rela- 
tively symmetric lags measurements extending to equally 
large velocities on both the red and blue sides of line 
— a clear indication of virialized gas motions, with the 
high- velocity line wings arising in gas closest to the cen- 
tral source (see R obinson & Perez 1990; Wels h fc Hornd 
ll991l : lPerez et al.l [l992: Ben tz et al.ll2009l for examples of 
how velocity resolved res ponses can be re lated to differ- 
ent BLR geometries; see iPetersonI ()200lD for a related 
tutorial). 

4. DISCUSSION 

The velocity-resolved reverberation signals exhibited 
by the H/3-emitting gas in the BLR of these three sources 
show striking differences in kinematic behavior. This is 
of particular interest because black hole masses derived 
from reverberation mapping are only valid under the as- 
sumption of gravitational domination of the BLR gas dy- 
namics by the black hole. The influence of gravity on the 
BLR in NGC 3516 and NGC 5548 is evidenced by the sig- 
natures of an infalling gas component in the former and 
a lack of significant radial motions of any kind in the 
latter (see Figure [3]). The apparent evidence for outflow 
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velocity (xlO^ km s"') 

Fig. 3. — Division of the H/3 rms spectral profile into equal-flux 
bins (top panels; vertical dashed lines), and corresponding velocity- 
resolved time-delay measurements (bottom panels) for NGC 3227 
(top), NGC 3516 (middle), and NGC 5548 (bottom), where the de- 
lays are plotted at the flux centroid of each velocity bin. Error 
bars on the lag measurements in the velocity direction (bottom 
panels) reflect the bin size, with each bin labeled by number in 
the top panels and negative velocities referring to blueshifts from 
the line center and positive velocities redshifts. Error bars on 
the lag measurements are determined similarly to those for the 
me an BLR la g (Peterson et al. 1998 with modiflcations described 
bv lPeterson et al.ll2004l V The horizontal solid and dotted lines in 
the bottom panel show the mean BLR lag and associated errors, 
as listed in Table 1, while the horizontal dotted-dashed line in the 
top panel represents the linearly-fit continuum level. Flux units 
are the same as in Fig. [T] 



in the case of NGC 3227 is a first for reverberation map- 
ping. In some sense, this should be no surprise given the 
overwhelming evidence for large-scale ma ss loss from the 
inner regions of AGNs (Crenshaw. Kra emer. fc Georg^ 
l2003f l. On the other hand, this does call into question 
the assumptions that allow us to estimate the mass of 
the central object, though we do hasten to point out 
that an outflow at escape velocity would still allow us to 
measure the black hole mass with the levels of accuracy 
currently claimed. Furthermore, the black hole mass we 
calculate for NGC 3227 of Mbh = i^-^tij) x IO^Mq, 
based on the lag in Table 1 and the line dispersion of 
the broad H/3 line measured from the rms spectrum (see 
D09b), is consistent within the statistical and system- 
atic uncertainties in this method to independent mea- 
surements using galactic stellar (Davics et al. 2006) and 
gas ()Hicks fc Malkanll2008i) dynamics. We also note that 
while the H/3 emission line is blueward asymmetric, the 
H/3 profiles of NGC 3516 and NGC 5548 are even more 
so, so we ascribe little importance to this. It strikes us 
as likely that we are observing a complex system such 
as a two-component BLR, in which there is an outflow- 
ing wind, in addition to a virialized, disk compo nent 
(|Murrav fc ChiandfijOTHEracleous fc Halpern|[2003l) . In 
this case, the mean lag that we measure, which is trac- 
ing the position of the majority of the line-emitting gas, 
arises from the disk component and is therefore under 
the gravitational influence of the black hole. Meanwhile, 
the high-velocity, blueshifted emission with very short 
lags arises from a wind at the inner BLR. Additional 
support for this scenario comes from the double-peaked 
profile shape of the rms spectrum, sugges ting a disk- 
like origin (e.g.. lEracleous fc Halperiilll994[ ). A further 
test for the virial nature of the BLR in NGC 3227 is to 
search for reverberation signals from multiple emission 
lines in this object. Data from this same campaign sug- 
gest that the HeiiA4686 emission-line flux also varied 
significantly over the course of the campaign, and future 
work is planned to search for a reverberation signal from 
this variable line emission. 

The observations of virial motions in NGC 5548 and 
the outflow in NGC 3227, in particular, are of fur- 
ther interest in the context of comparisons with previ- 
ous velocity-resolved studies of a siinilar nature (e.g., 
iSergeev et all 119991: [Poroshenko et al.ll2008f) . These are 
largely focused on NGC 5548 (the object for which we 
have the most reverberation mapping data), and a com- 
prehensi ve summary of many of the se past studies is 
given by iGaskell fc GoosmannI (|2008[ ) , who discuss the 
support for infalling BLR gas as implied by these vari- 
ous results, particularly for low-ionization lines like H/3. 
This is in contrast to what we clearly see in NGC 5548 
and NGC 3227. It is entirely possible that the velocity 
fields in these regions could change as a function of accre- 
tion rate, luminosity state, or over dynamical timescales. 
Interpretation of these new results is problematic and 
probably will remain so with only a handful of examples 
at single epochs. At this point, generalizing from these 
few sources is premature. 

5. SUMMARY 

In this work, we have presented three clear cases of 
differing velocity signatures of H/5-emitting gas from the 
BLR of three nearby AGNs, demonstrating the diversity 
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and probable complexity of the kinematic s in this re- 
gion. Our abili ty with this work and that of lBentz et alj 
see also lBentz et al.ll2008[ ) to recover statistically 
significant velocity-resolved reverberation responses for 
multiple objects is principally due to successful comple- 
tion of campaigns in which high quality, homogeneous 
observations were obtained over long durations (i.e., mul- 
tiples of the reverberation lags) with a sampling rate 
that was high compared to the relevant timescales be- 
ing investigated. Velocity-resolved reverberation map- 
ping studies such as described here are the next step 
toward producing a velocity-delay map, which will re- 
construct the two-dimensional kinematic structure of the 
BLR. This, in turn, will allow further insight into the ge- 
ometry and dynamics of the BLR, potentially leading to 
estimates of its inclination and ultimately reducing the 
systematic uncertainties in Mbh determinations. In par- 
ticular, placing direct observational constraints on the 
valu e of the reverberat ion mapping mass scale factor, / 
(see lOnken et al.][200l , even in individual objects, could 
reduce the scatter in the MRH-gy^ relation for AGNs 
jGebhardt et al."2000'; 'Ferrarese et~aI]l200H lOnken et all 
|2004; Nelson ct al. 2004). Despite the diversity in BLR 
kinematic signatures seen here, we do not see any sys- 
tematic trends in the location of these objects on the 
Mbh-ct* relation, suggesting that any effects the differ- 
ent velocity fields of these objects have on their black hole 
mass estimates are within the current scatter in this re- 
lationship. More detailed information about the velocity 
fields of these objects is needed to say anything further 
on the effect these kinematic differences have on their 
black hole mass measurements. Future work will focus 
on revealing more structure in BLR velocity fields as we 



work to create velocity-delay maps of the H/3 emission 
in NGC3227, NGC3516, and NGC5548 from the data 
presented here. 
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